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We consider the process of double photoionization (DPI) in a rare gas atom as a two-step process, namely (i) photoion- 
' ization in an inner shell followed by (ii) the emission of an Auger electron from an outer shell. An angular correlation 
function for the two emitted electrons is defined by analogy with the theory of angular correlation in nuclear physics. An 
(~j ' expression is obtained for this angular correlation function by a statistical method which makes use of the density and 
^ , efficiency operators. The latter takes care of the attenuation of the probability of detection of an electrons due to the 
J ■ geometrical properties of the detector. Theoretical values of the angular correlation function are obtained for DPI in xenon 
QO ' and these are shown to be in good agreement with the experimental results given by Kammerling and Schmidt jl^] . 
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The theory of angular correlation was originally formulated in nuclear physics. This was done for two possible kinds of 
phyical situation. 

(i) Angular correlation of successive nuclear radiations emitted by a radioactive nuclear species, e.g. a 7 — 7 cascade. 

(ii) Angular correlation of successive nuclear radiations resulting from the more general case of a nuclear scattering 
process. 

The earliest paper on this subject was written by Hamilton Jl| for the case of a 7 — 7 cascade. Later on, Gardner 
^ ■ wrote down the angular correlation function using the wave functions of the states occurring in a cascade decay. Racah 
— 1 ' |J obtained thereafter a simple expression for the angular correlation function in terms of the angle between the directions 
of emission of two successive radiations emitted by a nucleus. He did not worry about the history of this nucleus, and his 
theory held regardless of whether the emitted radiations consisted of bosons or fermions. Subsequently, Biedenharn and 
Rose Q extended Racah's work to give a general form of the angular correlation function for successive nuclear radiations. 
About the same time, Coester and Jauch || gave a theory in the more general context of nuclear reactions. Their paper 
' gives a derivation of the angular correlation function which brings out the statistical nature of the problem. Later, Devons 
, and Goldfarb wrote a detailed review of angular correlations up to that point. 

In our problem a rare gas atom absorbs a photon having specified properties i.e. energy, spin etc. As a result, the atom 
t-H ' emits a photo-electron giving rise to a singly charged ion. This ion now de-excites by emitting an Auger electron. This 
, yields a doubly charged residual ion. 
' Note two basic features of this problem. 

' (i) The initial atomic state is randomly oriented ( an improper state or "Gemisch"). 

(ii) Only limited information is available about the emitted electrons, usually their directions of motion. 
The density matrix or statistical matrix p was introduced into quantum mechanics to provide for the discussion of just 
\ such a situation |t], |l5| . It forms an important ingredient of the nuclear theory of angular correlations. It will naturally 
f> |- figure in our atomic theory of angular correlations as well. 

In section 2 we shall try to construct an angular correlation function for the atomic problem outlined above. This is 
. i-H , done in section 2.1 using a statistical approach . In section 3 we report our calculation for the xenon atom which was 
■ experimentally studied by Schmidt and his co-workersjl^j. 



2 An angular correlation theory for atoms : the double photoionization 
problem 

Double photoionization (DPI) occurs when an atomic target consisting of rare gas atoms is irradiated with photons from 
an advanced light source, e.g. a synchrotron, and an atom emits two electrons in quick succession. In case the time interval 
between the emission of these two electrons is very small, one could think of DPI as a one-step process. A considerable 
amount of work has already been done on this |[l4|l3| . However, in case the time interval between the successive emission of 
the two electrons is substantially longer than the time taken by the photo-electron to leave the interaction zone, DPI may 
be regarded as a two-step process ||ll|. This in its turn will depend on the energy imparted to the atomic system by the 
incident photon. 



1 



We shall now try to construct a theory for the two-step process described above, keeping in view the earlier work done 
in nuclear physics. Since this is a two-step process mediated by electromagnetic interaction at two vertices, we expect it 
to be a second order process as compared to normal photoionization (PI). Hence the probability of its occurrence will be 
much lower than that of PI, and its detection will call for much greater precision. 

Consider a randomly oriented rare gas atom in a S e state. It absorbs a photon and after a certain time interval emits 
a photo-electron from an inner shell, giving a single-vacancy ionic state. This intermediate ionic state now de-excites by 
emitting an Auger electron, typically from an outer shell, giving rise to a two- vacancy final state |0[, 
So our process amounts to 

hv + A — > A + + er — > A ++ + e± + . (1) 

The initial state of our system is characterized by the set of quantum numbers (J a M a a a ), or by virtual quantum numbers 
(J' a M' a a' a ) arising from possible interaction with other atoms and electrons. Here (J a M a ) or (J' a M' a ) are angular momentum 
quantum numbers, and a a ,a' a denote the set of remaining quantum numbers characterizing the system. 

We wish to obtain an angular correlation function for the two emitted electrons in terms of the angle 6 between their 
directions of emission. 



2.1 The method of efficiency and density operators 

We now approach the problem of angular correlation from a statistical standpoint. We note that there is a certain probability 
for the atom to attain the final state given above. This is described by the appropriate matrix element of the density or 
statistical operator p as defined in the literature |^,^]. Now, even if the atom goes over to the final state, because of the 
finite size of our detecting equipment and other limiting factors, this event may or may not be detected. There is thus a 
certain probability e (0 < e < 1) that the event will be detected. This probability is represented by the efficiency operator 
e which will depend on the size, position and geometrical arrangement of the detecting equipment. 

Now, how does the angular correlation function relate to the operators p and el We shall try to answer this question in 
the following manner. We start out by defining 

£ =E q cq(Q\p\Q) 

= J2q £ QPQQ ( 2 ) 
= Tr(ep) . 

Here eq is the efficiency or probability of detection of the state described by quantum numbers Q, and pqq the probability 
of the system being in the particular state Q. 

From the elements of statistical mechanics, we know that e is the expectation value ( or average value) of the efficiency 
operator e We shall presently see that the angular correlation function is simply related to e. 

Since e is the trace of a matrix it is invariant under a unitary transformation in Hilbert space. 

We now ask the question : what are the transformation properties of the matrices e and p as we go from one unitary 
representation to another? It is easily seen that e and p are tensor operators. Hence they are also called the efficiency and 
density (or statistical) tensors respectively [||. The fact that e and p are both tensor operators makes them amenable to 
further analysis. 

Since the initial atomic system is randomly oriented, we have rotational symmetry and angular momentum is conserved. 
Hence our state vectors are eigenvectors of J 2 and J z . In such a case the matrix elements of a tensor operator have a simple 
geometric dependence on the magnetic quantum numbers. This is given by the Wigner-Eckart theorem. 

Using this theorem we write the matrix element of the density operator as 

(J a M a a a \ p\f a M' a a' a ) = ^ (- 1 ) Ja_MaC M o M^ a ' 0fc « K »(" 7aQa ' J 'a a 'a), (3) 

where k a = J a + J a 'and K a is the projection of k a - 

Similarly, the matrix element of the efficiency operator is given by 

(J a M a a a \ e | J' a M' a a' a ) = 2J (-l) J "" M 'C^° a £i aKa (J a Q a , J' a a'a)- (4) 

k a K a 

Using the unitarity property of Clebsch-Gordan coefficients we get the density tensor of rank k a with (2k a + l) components 
pk aKa {J a a a , J' a a' a ) = 2_j (-l) Ma_M °C J J ^*° (iii (J a M a a a \ p\j' a M' a a' a ) . (5) 

M a M' a 

A similar expression can be obtained for the components of the efficiency tensor. 
So the expectation value of e becomes 

e = Tr(ep) = ^ p kaKa {J a a a , J' a a a )el aKa {J a a a , J' a Q.' a ). (6) 

J a J' a ct a a' a k a K a 

Our choice of phase ensures that pk a n a is a Hermitian tensor. It satisfies the relation 
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Pk aKa (JaOia, J' a a' a ) — (-l) Ja Ja + Ka Pk a - Ka (JLa'a, JaCHa)- (7) 

Some simplification yields the result 

pk aKa (Jaa a ,J' a a'a) = " i " ~ < J =II J2 ||Jb> II Jb)* {Jh\\ jl \\Ja) {M\ 3l II J °)* ■ (8) 

V 2Ja + 1 

The first decay, namely the photoionization process, is characterized by the relation 

Ja=Jb+jl. (9) 

Since the detection of the photo-electron and that of the singly charged ion are independent events, the joint probability of 
their detection is given by the product of the individual probabilities. Hence we can write 

{J a M a a a \ e \J' a M' a a' a ) = £ {J b M b a b \ e\J' b M' b a' b ) {jifii\e |jVi) 

vr .*J b JlJa f^Jbj'^L (10) 

Here (J b M b ), (J' b M' b ) describe the intermediate state and a b ,a' b denote the remaining quantum numbers. The summation 
in Eq. (0) extends over M b M' b n\[i' x . 
Then 

Sk aKa (J a a a , J'aa'a) = J2k 1K1 k bKb ek b^b( J i> a b, Jia' b )e klKl (jij[) 

il Ja } /-.-.n. 

«b«i*a./ri ~. att/ ~. ajx: rr. An: rr J r> r' I 

K fa K l K a 

with the relations 

J a = Jb + ji , k a = J a + J a ' , k b = J b + Jb' and ki = ji + ji'. 
Since the intermediate singly ionized state decays into a residual doubly ionized atom and an Auger electron, we can factorise 
the efficiency operator of the singly ionized atom in terms of the efficiency operators of the residual doubly ionized atom 
and the Auger electron. 

Since the residual doubly ionized atom is left in a sharp eigenstate with the quantum number J c and no further 
measurement is made on it, we get the efficiency operator 



J b 
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Ja 
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£k a K C (Jc<y c , J' c o;' c ) = v2JT+lSk c oS Kc oS Jc j' c S ai . a / a . (12) 
A matrix element of the efficiency operator for each electron has the form 



with i=l,2. i = 1 denotes the photo-electron and i = 2 the Auger electron. 

Using the axis of the detector as the quantization axis for each detected electron, we get 

(Qicn \jnn) = ^2 (° a ' \0i K i) Di&iWi)- (I 4 ) 

D(5Ri) is the corresponding rotation matrix for the ith electron. 
Then Eq. @ becomes 

SkiKiijij'i) =^ c k z K>.(jiji)D k K ** K ,('fti), (15) 

where 

ckiM'i)= 9FE WM (-i)^~ Mj <oo-ili*w>* 

x(0a' % \M)C^ iRi x (<T t \e\ai). (16) 

The symbol Ss represents a summation over the spin of the emitted electrons and depends on the characteristics of the 
detector. We now take the summation ^sf(aa') to be equivalent to E (a\ e \cr ) f(aa') . 

Making use of properties of Clebsch-Gordan coefficients and 9 — j symbols, and also using the additivity of rotation 
matrices p7| 

£^ K1 (3ii)^* 2 (3? 2 ) = Eu^nM^Q^ 1 ) n7 , 

we get 



3 



e = (2J b + l)(-) J ^- 2J ^(-l) k -^w(J b Jlj 1 j[,kJ a ) 

XtoCJb^'aia.WcJcfcMOiiiJ^CjaiapJ^^iflafe) (18) 

X (J c || ja lh) (Jc|| A \\J b }* (J b \\ jl \\J a ) (Ji|| ji || Ja>* • 

In Eq. ( |li| ) the summation is over ji, jl, ja, ja> &i «i andK2 . 

If the electrons are unpolarized, or if the detectors are insensitive to polarization, m = k 2 = and 000(618283) — 
Pk(cos82) . We now identify the angle 82 with 8, the angle between the directions of emission of the photo-electron and the 
Auger electron. 

From the Hermitian character of the efficiency and density matrices, it follows that 



ck.Uf) = (-i) K cL(i'i)- (19) 



For our case Ki = «a = 0, and 



We can thus write, 



CkoUf) = V¥+T^7+T ( _ ir -i + , c ^ o (20) 



g = (-i) J -+ J «- aj » v ^+TX) fc (-l)*- il - Ja «;(J6JJiiii,*J.) 

Xw(J b Jij 2 j2,kJo) {Jc\\j2 \\Jb) {Jc\\j 2 II Ji)* (21) 
X (Jill jl Pa) (Jill jl l|Jar c M(jlii)CfcoO'2j2)-Pfc(cOs6»). 

Here & is an even integer ranging from to k max , k max being defined as follows. Let {{ji + j[} max , {32 + 32} max } rnin = P- 
Then fc m<M = p if p is even and = p — 1 if p is odd. 

We now express (J c || J2 \\J b ) in terms of (J b \\j2 \\Jc) ■ Though the reduced matrix elements are neither real nor Hermitian, 
it happens that 

V2J„ + l(Ja|| j||J 6 ) = (-l) Ja - i+J W2J a + 1 (Ji|| 3 pa)* ■ (22) 

This gives 

e = E fe (-l) J1+J2c feo(jii 1 )^o(i2j 2 ) 

x (Jell 32 II Ji) (Je||i 2 II Ji)* (Jill jl II Ja) (Jill ji || Jo)* (23) 
xw(J b J b j 1 j[, kJa)w(J b J b j232, kJ c )Pk{cos8). 

If the finite size of the detector is taken into account, the efficiency of detection described by the matrix element of the 
efficiency operator must be changed slightly. Then we have to introduce Zk as the attenuation factor due to the finite size 
of the detector. We assume the detector to be axially symmetric Jic| , |l^ |. The efficiency tensor described by Eq.(fl5|) is now 
written as 

So the expectation value in our case becomes 

5 = E fc ^(l)^*( 2 )(-l) il+i2 Cfco(jlj04o(jaj2) 

x (Jc|| ja || Ji) (Jc|| ja II Ji)* (Jill ji II Ja) (Ji|| ji || Ja)* (25) 
Xw(J b J b jij[, kJ a )w(J b J b j2j 2 , kJ c )P k (cos8). 

Note that the 8 dependence of e is contained in the function 

w(0) =E k ^Wkz k (2)(-iy^c k0 (j 1 j' 1 )ct (j 2 3 , 2 ) 

Xw(J b J b jij[ , fcJa)w(J 6 Jjj 2 j2, k.J c )P k {cOs8). 

We now define W(8) to be the angular correlation function for the two emitted electrons where 8 is the angular separation 
between their directions of emission j9[. It is clear that angular correlation between the directions of emission of the photo- 
electron and the Auger electron is a direct manifestation of the efficiency of the observing equipment. In section 3 we shall see 
that the angular correlation function so defined agrees closely with the measured angular correlation in the DPI experiments 
on xenon jljj]. This confirms that the observing equipment does play a role, introducing an element of probability which 
finds expression in the angular correlation function. 

This definition will have to be modified if it is possible for the photo-electron to be emitted into more than one angular 
momentum channels. In section 3 we shall see how this modification can be made. 

Note that Eq.(|2^) holds formally not only for double photoionization in atoms, but generally for two-step angular 
correlation experiments involving either fermions or bosons. 
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3 Calculation and results 



We consider the problem of DPI in xenon. Neutral xenon atoms are irradiated with a photon beam of energy 94.5 eV. This 
leads to photo-ionization in the 4d 5 / 2 shell followed by an iVs — 02,302,3 So Auger decay. Using the dipole approximation 
the possible photoionization channels are e)4d 5 / 2 — ► ep/7/2 , f)4ds/9 — > ep/5/2 and g)id 5 / 2 — ► £j>f>3/2 respectively. The 
Auger transition is characterized by only one partial wave £Arf5/2 |18| . These transitions are governed by the corresponding 
selection rules for photoionization and Auger transitions. 

Since the initial photoionization process is not characterized by a single angular momentum quantum number but by 
angular momentum quantum numbers corresponding to three possible channels, the angular correlation function described 
in section 2.1 above will be modified. The total intensity will, however, remain unchanged. If the photoionization channels 
are described by the total angular momentum quantum numbers jf, j/ and jf and the Auger electron by j 2 , then the 
expectation value of the efficiency operator becomes 



e(9) = e.(0) + E t {6) + e g (9) + e ef (9) + e fg (9) + e g£ (9) . (27) 

Here 



X \{Jc\\j2\\J b )\ 2 \{Ja\\jt\\J b )\ 2 (28) 

xw(J b J b jfjl, kJ a )w(J b J b j2j2,kJ c )P k (cos9). 

The expectation values £/(#) and £ g {9) have the same form with jf — ► j{ and jf — > jf respectively. The quantity £ij(6) 
is an interference term arising from interaction between photoionization channels i and j = e, /, g with i j) . 



£e/(0) = E fc ^(l)2fe(2)(-l) J2 |(Jc|| 32 \\Jb}\ 2 <Ja|| \\J b ) (Ja\\j( \\J b ) 

x[(-l) jt c k0 (jtj[) + (~iy'c k0 (j{jt)]ct Q (j2j2) (29) 
xw(J b J b jfj{,kJ a )w(J b J b j2j2, kJ c )P k (cos9), 



= E k ^(l)z k (2)(-y* \{J c \\j2\\J b }\ 2 (Ja\\j{ ||J»>(J,||jf ||Ji> 

x [(-l^cfcoO'i'if ) + (-i) il9 c fco (j?jf)]4oOM2) (30) 

xw(J b J b jfj{, kJ a )w(J b J b j2j2,kJ c )P k (cos9) 



and 



= E fc %(1)«*(2)(-1)» |(Jc||j2 ||Jb>r \\Jb) (Ja\\fl \\Jb) 

x [(-l)^o(jfjf) + (-1)'? «%o(j'f Ji)]cJo(jaia) (31) 
xw(J b J b jfjf, kJ a )w(J b J b j2j2, k,J c )P k (cos9). 

In £ e (#) ,£f(9) and £ 9 (0) fc is the smallest even integer of the sets {2jf , 2j 2 , 2J b }, {2j{,2j 2 ,2J b } and {2jf , 2j 2 , 2J b } respec- 
tively. In s e f(6), £f g (9) and e ge (9) k is the smallest even integer (k 7^ 0) of the sets {jf + j{ , 2j 2 , 2J b \, {jf + jf , 2j2, 2 J 6 } 
and {jf + jf ,2j2,2J b } respectively. 
We can now write 



where 



Here 



£W = \(Jc\\j2\\J b }\ 2 [iJe(9)+5 2 LU J {9)+S 2 2 LU g (9)+8 1 U ef {9)+S2U Ja {9)+S 1 82LU ge (9)] , (32) 

x (Ja\\3i\\ J b) „„j X (Mill \\Jb) 

dl - I T II -c II 7 \ ffln » d 2 = / 7 II - e II T \ • (33) 
Ja J Ji, (Ja\\3t\\Jb) 



Ve(0) = E k Ml)^)(-lY 1+32 Ck (jtjt)cl (j2j2) 

o{JbJbji3x, kJ a )w(J b J b j 2 j2,kJ c )P k (cos9) 



xwl 



(34) 



LOf{9) and u> g (9) have the same form with jf — > j/ and jf — > jf respectively. 
And 



co ef (9) = z k (i) Zk (2)(-iy-[(-iytc k0 (jtj() + (~iy(c k0 (j(f 1 )]ct (j2j2) 

Xw(J b J b jfj{, kJ a )w(J b J b j2j2, k.J c )P k (cOs9), 

= E**k(i)%(2)(-i)^[(-i)^c w ,0'i / if) + (-i) i? c fc o0Tjf)]4 (j 2 j 2 ) 

xw(J b J b jfj{ , kJ a )w(J b J b j 2 j2, kJ c )P k (cos9) 



(35) 
(36) 



and 



= E^ai)^42)(-iyn(-l)^C fc o(jf^) + (-l)^ Cfc o(jfjf)]4o(j2j2) (37) 

xw(J b J b jfjf,kJ a )w(J b J b j 2 j2, kJ c )P k (cos9). 
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The parameters <5i and S 2 can be determined uniquely by comparison with experiment |(|. 
Writing out the expression (|32|) in terms of Legendre polynomials, we get 



e{6) = N [a + a 2 P 2 (cos9) + (IaPa{cos9)\ . (38) 

Here 

N=|<J c ||j2||J i) )| 2 |(J £l ||j 1 e ||J b )| 2 , 

and the coefficients are a = z (l)zo(2)(l + <5? + <5 2 ) , a 2 = z 2 (l)z 2 (2)(1.0204 + 0.751<5j ! +0.852 + 0.269<5 2 +0.311<5 a + 0.733<5i<5 2 ) 
and a 4 = z 4 (l)z4(2)(0.5510 - 0.1225?). 
This gives, 

s = W [1 + b 2 P 2 (cos9) + b A Pi{cos9)] , (39) 

where 6 2 = = 0.760, 6 4 = = 0.042 and K' = Ka = |(J C || j 2 || J 6 )| 2 K^Hi? ||Jb}| 2 o . 
We now define the angular correlation function for this case of channel mixing by writing 

W{9) = [1 + b 2 P 2 (cos9) + b i P i (cos9)] . (40) 

Figure 1 gives the results of comparison between our theoretical values and the experimental values It will be seen that 
our polar plot of W{9) agrees quite closely with that given by Kammerling and Schmidt jl4| except for a small difference 
in the region around 9 = degree and 9 = 180 degree. 

By using scattering theory with the appropriate boundary conditions, it is possible to obtain <5i and 5 2 without recourse 
to the experimental curves. This in its turn should give not only oo, but also b 2 and 64. This is the way our theory can be 
used to predict the value of the angular correlation function at any required angle. This would make our theory autonomous. 
In the absence of multichannel interaction, 5\ = 82 = 0, Eq.dlc]) reduces to Eq.(p6|). 

The expectation value e turns out to be the product of a normalization factor depending on the reduced matrix elements 
and an angular factor. The simplicity of the latter is very striking. It seems as if the dynamical calculation involving the 
radial matrix elements is redundant. Of course this is not true, because to extract dynamical properties of the system such 
as triply differential cross sections (TDCS) we have to play with the normalizing factor. However, the very simplicity of the 
result hinges on the factorization of the problem into a dynamical part and an geometrical part depending on 9. This comes 
from the use of Wigner-Eckart theorem, which is a consequence of the fact that we are dealing with the matrix elements of 
tensor operators. 
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4 Appendix A 

The following is a description of the notations used in this paper. 
A : neutral rare gas atom which undergoes double photoionization. 
e,\ : photo-electron. 
: Auger electron. 

J a , Jb, J c : total angular momentum vectors of initial, intermediate and final states of the rare gas atom respectively. 

J a , Jb, Jc '■ total angular momentum quantum numbers of initial, intermediate and final states of the rare gas atom respec- 
tively. 

M a , Mb, M c : projection quantum numbers for the initial, intermediate and final states of the rare gas atoms respectively. 

oia, otb, ce c '■ other quantum numbers corresponding to initial, intermediate and final state of the rare gas atoms respectively. 

ji, J2 : total angular momentum vectors of the photo- and Auger electrons respectively. 

31 , j 2 : total angular momentum quantum numbers of the photo- and Auger electrons respectively. 

/J.i,fl 2 '■ projection quantum numbers corresponding to ji and j 2 respectively. 

h,l 2 : orbital angular momentum quantum numbers of the photo- and Auger electrons respectively. 
s\,s 2 : spin quantum numbers of the photo- and Auger electrons respectively. 
<7i,<7 2 : spin projections of the photo- and Auger electrons respectively. 
p : density operator. 

pk iK ,i '■ density or statistical tensor of rank fcj with (2ki + 1) components, m being the projection quantum number denoting 
one component. 

e : efficiency operator. 
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EkiKi ■ efficiency tensor of rank ki with (2fcj + 1) components, m being the projection quantum number denoting one 
component. 

e : probability of detection of any event. 

3 \\Jf) '■ a typical reduced matrix element. 

s : expectation value of the efficiency operator. 

Qi : angular direction of the i th electron. 

w : Racah coefficients. 

Cm{m 2 m '■ Clebsch-Gordan Coefficient obeying the triangle rule Afjijzj) and mi + m,2 = m. 

9 : summation over all non-observed properties of the atom and the electrons, e.g. spin states and characteristics of the 
detecting equipment. 

Ck K : Attenuation factors corresponding to the state of polarization. 

Zfc(l) : attenuation factor due to the finite size of the detector which detects the photo-electron. 
Zfc(2) : attenuation factor due to the finite size of the detector which detects the Auger electron. 
e p : energy of the photo-electron. 
£A : energy of the Auger electron. 

All primed quantum numbers denote virtual states. 
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Fig 1: A polar plot of our angular correlation function for DPI in Xenon is compared with the experimental polar plot 
given by Kammerling and Schmidt |14J. 
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